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Abstract 

It is shown that the logarithmic derivative of the characteristic polynomial of a Wilson 
loop in two dimensional pure Yang Mills theory with gauge group SU(N) exactly satisfies 
Burgers' equation, with viscosity given by 1/(2N). The Wilson loop does not intersect 
itself and Euclidean space-time is assumed flat and infinite. This result provides a precise 
framework in 2D YM for recent observations of Blaizot and Nowak and was inspired by 
their work. 
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1 Introduction. 

Recent numerical work provides evidence that Wilson loops in SU(N) gauge theory in 
two, three and four dimensions exhibit an infinite N phase transition as they are dilated 
from a small size to a large one; in the course of this dilation the eigenvalue distribution 
of the untraced Wilson loop unitary matrix expands from a small arc on the unit circle 
to encompassing the entire unit circle [TJ [2] . An analogous effect takes place in the two 
dimensional principal chiral model for SU(N) [3]. 

The universality class of this transition is that of a random multiplicative ensemble 
of unitary matrices. The transition was discovered by Durhuus and Olesen [2J (DO) 
when they solved the Makeenko-Migdal [5] loop equations in two dimensional planar 
QCD. The associated multiplicative random matrix ensemble [6] can be axiomatized 
in the language of noncommutative probability [7]. It provides a generalization of the 
familiar law of large numbers. The essential feature making a difference is that one 
case is commutative and the other not. Various recent insights into the DO transition 
El [TO] point to possibly deeper interpretations of the transition. 

In this note, motivated by a recent paper by Blaizot and Nowak [10J, I present an 
exact map from the average characteristic polynomial associated with a Wilson loop to 
Burgers' equation. This extends to finite N the original work of DO at N = oo, where 
the inviscid Burgers' equation plays a central role. The main observation is that all finite 
./V effects are exactly represented by reinstating a finite viscosity in Burgers' equation, 
given by j^. Positive N gives positive viscosity, so the equation knows at least that N 
should not be negative. I suspect that integral iV's are identified as special by a Mittag- 
Leffler representation of the solution, stemming from a product representation of 
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the average characteristic polynomial, and depending also on the initial condition. 

In addition to making the insight of [10J particularly transparent, I hope that this 
result would also aid future efforts to exploit large N universality in dimensions higher 
than two for obtaining analytical quantitative estimates of the ratio between a scale 
describing perturbative phenomena and the scale of confinement. This was the original 
motivation for seeking to establish numerically large A?" phase transitions in Wilson 
loops pQ. 



2 Characteristic polynomial. 

An N x N simple unitary Wilson loop matrix W, defined on a curve that does not self 
intersect, with r denoting the dimensionless area in units of the 't Hooft gauge coupling, 
has the following probability distribution: 

P N (W, r)dW = £ d RXr (W)e- TC ^dW (1) 
R 

The sum is over all irreducible representations R with character Xr(W) anc ^ second order 
Casimir C2{R)- dW is the Haar measure. Normalization conventions are standard [2] 
and r > 0. We introduce the average characteristic polynomial 

Q N (z,r) = (det(z-W)) PN{T) (2) 

One can think about Qn( z , t ) as the generating function for the (xr(W)) with totally 
antisymmetric R. Simple manipulations [2] produce an integral representation: 
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Q n (z,t) = [°° due-%™ 2 \ z - e- T( - u+1 M] N (•'>) 

It is more convenient to study 

q N (y,r) = (-l) N Q N (- e y , r) (A) 



where, for the time being, y is kept real. qN(y, t) is even in y and this is the main reason 
for extracting the exponential factor from Qn- Changing the integration variable u to 
x = y + t(u + 1/2) gives: 



qN (y ]T) = JJL I ~ dxe -£(y- X )\Nlo g (2cos Hx /2)) 



3 Main result. 

It is now a trivial matter to observe that 

dqN _ 1 d 2 q N 
dr ~ 2N dy 2 



(6) 
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with initial condition 

lim[ 9jv (y,r)] = (2cosh(y/2)) iV (7) 

The behavior at y — > ±00 prevents solving ([6]) by Fourier decomposition and any asso- 
ciated general conclusions about boundedness as r — > +00. The initial condition is a 
consequence of 

Pjv(W,0) = 5{W,1) with y dW r (5(H / ',Wo)/(l^) = /(W ) (8) 

for any Wo £ SU(N). This equation can be also directly derived from the polynomial 
formula of Qjv> without going to the integral representation. This heat equation is 
related to Burgers' equation (for example, see [12], problem 12(a), p. 214) by 

1 / \ 1 9 log gjy(y,r) 

<M^) = -- ¥y (9) 

Burgers' equation and the initial condition are 

+ < PN^- = 7777^3-' 4>N{y,0) = -- tanh - (10) 



0r 8y 2N dy 2 ' rJVVy ' ' 2 2 

At N = 00, N drops out of the equation giving the inviscid limit: 

86 86 , , 

The initial condition is N independent so we can drop the N subscript on 6 at N = 00. 
So long as 6 is uniquely defined, this is the point- wise N = 00 limit of 6^. 

The equation can be solved by the method of characteristics (for example, see |12j . 
p. 16.) for an arbitrary initial condition 

6(y,0)=h(y) (12) 

The solution is given implicitly by 

6(y,r) = h(y -T6(y,r)) (13) 

This equation is known to produce a shock at a time r* > which is the first time at 
which multiple solutions become available, r* is the smallest positive value satisfying 

with (d 2 h/dy 2 ){y*) = (14) 



(dh/dy)(y*) 

We are interested only in solutions odd in y; hence, assuming h(y) to be smooth near 
y = we expand: 

h(y) = ay + by 3 + cy 5 + .... (15) 
This implies that y* = and therefore 

t* = -- (16) 
a 



A shock will form if a < 0. In the case of N = oo 2D YM we have 

Hv) = ~\ ^nh | = -y/4 + y 3 /48 - .... (17) 
Therefore, the critical area corresponds to 

r* = 4, (18) 

the well known critical value [U [6J . 

Universality can be invoked now in a sense that applies to the nonlinear equation 

producing a generic shock [13|. I14j . This means taking the simplest polynomial h(y) 
capable of producing shocks: 

h(y) = ay + by 3 (19) 

with a < 0, b > 0. The y location of the shock is at the origin, y = y* = 0. Extending 
h and y to the complex plane provides a geometric view of this universality in terms of 
the structure of the evolving Riemann surface y((ft, r) parameterized by r > 0. One can 
also take r into the complex plane. 



4 Large N = small viscosity. 

Making the viscosity nonzero is a singular perturbation which eliminates the shock and 
has the same effect as making N finite. Large N universality will hold in the vicinity of 
the critical area and corresponds to universal behavior in the vicinity of the would-be 
shock for small viscosities, which is the simplest dissipativ^] regularization of the shock. 

The important new insight is that the large N transition is equivalent to a mov- 
able singularity, determined by the initial condition, rather than by the evolution rulell 
Thus, the simplest initial condition producing a shock will also lead to a universal small 
viscosity smoothing of the shock. 

Running the derivation backwards, with the minimal initial condition 

h{y) = -y/4 + y 3 /48 (20) 

produces an integral representation on which a double scaling limit can be taken directly, 
exactly reproducing the limit used in matching to the large N transitions in higher 
dimensions than two in [H[2]. The critical exponents [i = 1/2,3/4 associated with the 
scalings iV 7 ' that need to be taken [2] are identical to those found in defining the small 
viscosity limit [IB] . The associated integral, studied in detail in [2] (see [T7] for a plot), 
is related to Pearcey's integral by a contour change, as indicated in |10j . 

The particular initial condition (|19p has been analyzed in great detail in |16j . 

1 The shock can be regulated also dispersively, in which case we could use a third derivative on the right 
hand side of the inviscid Burgers' equation, producing the KdV equation. If there were a symmetry restricting 
to a Hamiltonian partial differential equations, this might have been the equation defining the universality 
class. 

2 Somcthing similar happens in the context of models consisting of one or several large matrices, where 
Painlcve equations enter (see for example [5] and [T5]). 
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5 Higher critical points. 



We have become accustomed to expect higher critical points, of reduced degrees of 
stability, to accompany a basic large N critical point. Looking at (|15p it seems plausible 
that setting 6 = and making c > would produce a critical point of one degree of 
stability less. Obviously, if this is true, a whole hierarchy will be generated, by initial 
conditions of the form ay + by 2 " 1-1 with integer m > 2. If one is not worried about 
the convergence of the associated universal integrals and one is also willing to give up 
the y — > —y parity symmetry, also higher critical points with half integer m could be 
studied, at least as formal originators of asymptotic series. 

It would be intriguing if parent models existed with physical symmetries that selected 
one of these higher critical points. More work on this is left for the future. 



6 Product representation. 

It certainly is true that 

N 

Q N (z,t) = (det(z - W)) Pn(t) = l[(z - Zi (r)) (21) 

l 

One may view the Zi{r) as certain averages of the eigenvalues of W, but not as usually 
defined: 

N 

detC* -W) = ]J(z - Zi(W)), zf{T) = (zi(W))p N{T) (22) 
l 

In [2] it was proved that |^(t)| = 1 for all i = 1, ...,N (see [17] for a plot); this indicates 
that the Zi(r) are to be viewed as an approximations to the zf v (r). By applying large 
iV factorization N times, one can argue, at least away from large iV critical points, that 
identically ordered z£ v (tYs and Zj(r)'s are equal to each other. I suspect that this stays 
true also in the double scaling limit. If this suspicion is validated, we shall obtain a new 
method to identify, using numerical simulations, the location and nature of the large N 
transition in dimensions higher than two. 

It is therefore interesting to derive evolution equations for the z%(j\ After inserting 
the product (f2Tj) into the heat equation © and applying @), standard manipulations of 
the kind employed in the study of Calogero systems produce 

z 1= lyz_ h + z 1 

Here zj = dzj(T)/dr and J2k means that the k = j term is dropped from the sum, 
where the index k runs from 1 to N. This equation is form invariant under Zj — > 1/zj 
and expected from the structure of the polynomial. In addition, again as 

expected, the product of all zeros is constant in r. Moreover, the equations of motion 
m imply d\ Zj (T)\ 2 /dr = 0,j = l,..,N. 
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The initial condition is %(0) = 1, j = 1,..,N and is degenerate. However, at any r > 
the degeneracy is lifted; for example, at infinite t, we have Zj (oo) = e 2m (j+ 1 / 2 - N / 2 )/ N . 

The map z = — e y creates an infinite number of copies of the zeros Zj, which are all 
on the imaginary axis. We choose one specific yj for each zj,j = 1,..,N. The equation 
of motion for the j/j's is: 

y, - ±Z = ±£ E ^ < 24 » 

The universal description changes the equation obeyed by the y/s. However, as pointed 
out in [16] on the basis of an old theorem [18], the yj(r) still stay on the imaginary 
axis for all r. In the universal case periodicity under yj(r) — > yj(r) + 2mjiri, rrij £ Z 
is lost, since the initial condition on the y^s no longer is periodic. Thus, one needs to 
use the yj variables to make the connection between the exact equations of motion and 
the universal ones. I leave a more detailed study of the universal limit of the eigenvalue 
motion to the future. 



7 Large r behavior. 

The regularization of the shock provides a smooth connection between small and large 
loops. In two dimensions Burgers' equation provides an exact renormalization group 
type of equation allowing the evaluation of 4>n(v, t) when r — > oo, given <p^(y,T) in 
the limit r — ► 0. The approach to the limit r — > oo gives the dimensionless string 
tension associated with the dimensionless area r. Here we only show how the correct 
4>]y(y,T = oo) is obtained. It is clear that Qn(z,t = oo) = z N + (— 1)^. This simply 
says that at infinite r all {W m ) terms, for any m > 0, can be replaced by zero. 
Using (jlj), we conclude that the large r behavior is given by: 

lim ( e - — q N (y,T) )= 2cosh— U - (25) 

We now wish to recover the ensuing <p^{y,T = oo) from Burgers' equation. The 
route is again in reverse of our derivation: First go to the heat equation, then get the 
integral representation in order to incorporate the initial condition. Finally, in order to 
get the asymptotic behavior for large r, change variables in the integral representation, 
arriving at: 

jjdy logq N (y,T) = 

f due-^ sinh(( Mv /7+y)/2)(2 co S h((«V7+j / )/2)) Ar - 1 , . 

( 26 ) 

J due 2-(2cosh((M0 I +y)/2)) JV 

For large r, one of the two exponents making up each hyperbolic function dominates, 
depending on the sign of u: 



T li i"So(^^ 10g ^^ T )) 
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Here, e(u) is the sign function. The above equation implies that 



</W(y, oo) = lim 




) y logq N (y,T) 



itanh^ 



(28) 



This is the expected result. 

At infinite N, the hyperbolic tangent becomes a sign function. In an electrostatic 
picture it is obvious that the above result holds if the poles of (p(y, r) are uniformly 
spaced and dense on the circle \e y \ = 1: Viewing the poles as charges, the jump e(y) 
comes from crossing the line charge at z = —1 as y goes through zero along the real 
axis [2]. That the solution has this limiting behavior is essential for confinement, which 
would be indicated by the leading correction to the above result being exponentially 
small in r. 

Note that r was taken to infinity at finite N; the final result admits a subsequent 
infinite N limit. Had we taken N — > oo first, we could have interpreted the shock, 
appearing first at r = 4, as a jump between two extremal solutions of the implicit 
equation defining the solution for r < 4. With the wrong initial conditions this jump 
might not grow to the full size required for consistency with confinement; thus, the 
transition in itself is insufficient to guarantee confinement. If we want to add the input 
that there is confinement we need to put a constraint on the initial condition. 

Regarding [10], following [22], I opt not to address here the question how Burgers' 
equation relates to turbulence. As a start, I refer the reader to [23]. In general, one would 
hope that the analogy to the three dimensional incompressible Navier-Stokes equation 
does not hold too literally. Large N would map to large Reynolds numbers, while small 
N to small Reynolds numbers; however, I am hoping that matters simplify at large N - 
if they do not, one would be better off concentrating on N = 3. 

Again, I leave details for further work. 



The primary objective of this paper was the derivation of (|10p as an exact equation 
holding in two dimensional Yang Mills theory with gauge group SU(N) defined on the 
infinite Euclidean plane. A surprising simplicity in the area dependence of the average 
characteristic polynomial of simple Wilson loops was found. Nevertheless, the essential 
feature of the existence of a large iV phase transition is captured by this observable. In 
this respect the average characteristic polynomial of the Wilson loop is superior to traces 
of the Wilson loop in some fixed representation. As explained in [2] this observable has 
other advantages, in dimensions three and four. 

The simple and exact finite N relation to Burgers' equation presented above seems 
to provide opportunity for progress in different directions, as emphasized in the course 
of this paper. The secondary objective of the paper was to present enough observations 



8 Discussion. 
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to convince the reader that there are many interesting issues left to explore. Last, but 
not least, the insights of Blaizot and Nowak [10] deserve further study. 

The shock at r = 4 is reminiscent of the possibility that instantons at infinite N 
might herald, as r — > 4", a jump in certain particularly sensitive quantities in 4D YM 

EE]- 

It should also be mentioned that workers in lattice field theory [20] have shown 
numerically that in four dimensions the trace 2cos# of a Wilson loop for SU(2) seems 
to evolve with the area as if 9 were diffusing on the SU(2) group manifold where the 
eigenvalues of W are e =t * e . For N = 2 there is no essential distinction between the 
characteristic polynomial and any other gauge invariant observable related to the matrix 
W. 

9 Added note. 

Blaizot and Nowak [21] have independently identified the viscosity as 
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